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ABSTRACT 

Refraction and diffraction of incoming radio waves by the ionosphere induce time 
variability in the angular positions, peak amplitudes and shapes of radio sources, po¬ 
tentially complicating the automated cross-matching and identification of transient 
and variable radio sources. In this work, we empirically assess the effects of the iono¬ 
sphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We 
directly examine 51 hours of data observed over 10 nights under quiet geomagnetic 
conditions (global storm index Kp < 2), analysing the behaviour of short-timescale 
angular position and peak flux density variations of around ten thousand unresolved 
sources. We find that while much of the variation in angular position can be attributed 
to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz 
are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching 
under typical conditions. By examining bulk trends in amplitude variability, we place 
upper limits on the modulation index associated with ionospheric scintillation of 1-3% 
for the various nights. For sources fainter than ^1 Jy, this variation is below the image 
noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is 
not a significant impediment to the goals of time-domain science with the MWA at 
154 MHz. 

Key words: atmospheric effects — instrumentation: interferometers — radio con¬ 
tinuum: general 


1 INTRODUCTION 


The fast survey capabilities of next-generation radio tele- 
scopes such a s the Murchison Wid e field Array (MWA; 
Lor^dale et al.l l2009l : iBowman et al.l l2013l : iTingav et alJ 


2013ll , th e Australian Square Kilometre Array P athfinder 
(ASKAP; I Johnston et al.l [20081 : iHotan et al.l 201^}and the 


Low Frequency Array (LOFAR; Ivan Haarlem et ^ l2013l l 
greatly facilitate the study of transient and variable radio 
sources. The primary observables in such surveys are sky po¬ 
sitions and flux densities, while secondary observables may 
include angular source sizes and shapes, polarisation frac¬ 
tions, polarisation vectors and rotation measures. 

Many astrophysical phenomena, including gamma-ray 
bursts, stellar flares, pulsars, exoplanets, black hole jet 
launches, supernovae and fast radio bursts produce tran¬ 
sient and variable radio emission on a wide range of time 
scales, from milliseconds to days. Locating and identifying 
the progenitors of these events has been the focus of many 
recent blind surveys making use of new low - frequency tele¬ 
scope s lljaeeer et al.l [ioi^: iBell et al.l l20l4 ICarbone et al.l 
I2OI4I : lObenberger et a l.||2015|) and also upcoming surveys at 
higher frequencies ( Murphy et al.ll2013lL Given the high data 
rate from these new instruments ( Norrisir201(jl . the extrac¬ 
tion, cross-matching, light-curve generation and classifica- 
tion of sources must be conducted in an automated fashio n 
llMurphv et al.|[2013l : ILo et al1l20l4 ISwinbank et al.ll2015h . 
An important aspect of these automated pipelines is that 
major sources of error and extrinsic/foreground variability 
are understood and accounted for to prevent source mis- 
classificat ions. The recent discove ry of radio emission from 
fireballs llObenberger et ^ |2014^ exemplifies a previously 
unknown terrestrial source of transient radio emission that 
may confuse searches for celestial transients. 

One important contributor to extrinsic variability, es¬ 
pecially at frequencies less than ~lGHz, is the distortion 
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of phase fronts due to propagation through the Earth’s 
ionosphere. The ionosphere is the ionised component of the 
Earth’s atmosphere. It extends between altitudes of roughly 
50-1000 km, with the ele ctron density peaking near 300 - 
400 km altitude at night l|Luhmannl 1 19951 : ISolomonI I2OI0I L 
At low radio frequencies, spatially-varying electron den¬ 
sities in the ionosphere induce variations in the angular 
position and flux density of radio sources throu gh refrac¬ 
tion and diffraction of propagating radio waves llBougeren 
Il98ll : Ijacobson fc Ericksonl1l992bl : IXassim et ^l2007ll . Po- 

sition shifts can diminish the accuracy of cross-matching 
between different epochs of observation, and propagation- 
induced amplitude variability can result in celestial source 
misclassifications. Smearing and angular broadening can 
cause otherwise point-like sources to become distorted, and 
this can pose difficulties for deconvolution and source iden¬ 
tification/extraction. Quantifying these effects is important 
both for the design and optimisation of algori thms that re¬ 
move ionospheric distort ions from the data (ICotton et al.l 
I2OO4I : lintema et al.l l2009h , and also to obtain realistic er- 
ror bars on ob servables for time-domain astrophysics (cf. 
iBell et al.ll2014h . 

A number of studies have been aimed at quan¬ 
tifying the effects of the ionosphere on radio astro¬ 
nomical observations at low radi o freque n cies. These 
include the t heore tical work of iHewis 3 (1195111 and 
iMever-Verned lll980ll . and observational work with th e 
Nangay radio interferometer llBougeredll98ll : lMercie ^ra¬ 
the Westerbork Radio S ynthe s is Telescope dSpoelstral 
1 19831 : ISpoelstra fc Keld^ Il984l : ISpoelstr^ Il985l. 1997), 
the Los Alamos radio interferometer ll Jacobson et al.l 


Array 

WLA) radio interferometer djacobson & Erickson 

1992allb 

, I993I: Hoogeveen & Jacobson! 1997al: Coker et al. 

2OO9I: ICohen & Rottgerins 

2OO9I: Helmboldt et al. 2012a bl: 

Helmboldt & Intemall2012l:lHelmboldt||2014|). These have es- 


tablished that radio telescopes are highly sensitive to and 
can even be used to characterise ionospheric perturbations 
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and irregularities ordinarily studied through other means, 
by isolating propagation-induced fluctuations in the ampli¬ 
tudes, phases and angular positions of radio sources. The ef¬ 
fects of electron density fluctuations occurring on timescales 
of several minutes and with spatial scales of ^lOkrn have 
been observed in interferometric data. This is important be¬ 
cause the characteristic spatiotemporal scales of ionospheric 
behaviour determine the spatiotemporal scales on which cal¬ 
ibration is required, and ultimately the pr ecision to wh ich 
astronomical observables can be measured llBraunll2013h . 

The types and occurrence rates of ionospheric fluctua¬ 
tions depend on geographic location, tropospheric weather 
and geo magnetic condition s (ISchunk fc Nagyll2009|l. Recent 
work by IPatta et al.l (I20l4l and Arora et al.l ( 2015ll has ex¬ 
plored the possibility of using Global Positioning System 
(GPS) satellite measurements to calibrate for ionospheric 
effects in radio interferometric data. However, the exis¬ 
tence of a substantial population of irregularities on length 
scales below the resolut ion of global GPS maps (~100km; 
iRideout fc Cost^ I 2 OO 6 I I implies that GPS measurements 
alone cannot fully predict the extent to which the iono¬ 
sphere affects radio observations dLoi et al.ll2015al l. Further¬ 
more, the coherence properties of ionospheric irregularities 
of a given scale size depend on the configuration of the 
array, with decoherence and scintillation effects more se¬ 
vere on longer base lines llLonsdalell2005l : [Xassim et al.ll2007l : 
lintema et al.l 1^091 1. Although general statements can be 
made based on physical arguments and established knowl¬ 
edge (see lj2]), a direct examination of the data is the most 
straightforward means of assessing the severity of iono¬ 
spheric effects for a given radio interferometer and geograph¬ 
ical site. 

This work aims to empirically establish the effects of the 
ionosphere on radio astronomy conducted with the MWA, 
focusing on its impacts on time-domain science, for which 
the two most important observables are angular position and 
flux density. We directly examine the statistics of those ob¬ 
servables in 10 nights of data collected as part of the MWA 
Transients Survey (MWATS; PI: M. E. Bell), a general- 
purpose, low-frequency (154 MHz) radio survey looking for 
transient and variable sources. We begin in )J2] with an 
overview of how the ionosphere affects radio interferometric 
measurements. We describe our observations and analysis 
approach in SjS] and present the results in Sjd] We interpret 
the results and discuss their implications in an 

Angular broadening and shape distortion effects are not 
considered here, although they may diminish the perfor¬ 
mance of astronomical source-finding algorithms (by making 
the source more difficult to identify, or increasing the mea¬ 
surement error in flux density and/or position). As argued 
in i)2.1l distortion and broadening effects are expected to be 
secondary to angular position shifts of whole sources, given 
the compact physical size of the MWA, for sufficiently short 
integrations. 


2 OVERVIEW OF IONOSPHERIC EFFECTS 

Electron densities in the ionosphere are commonly described 
in terms of the total electron content (TEC; i.e. the electron 
column density J Al, where rie is the local electron num¬ 
ber density and I parametrises distance along the line of 


sight), whose units are electronsm“^, commonly expressed 
in TEC units (TECU; 1 TECU = 10^® electrons m“^). Radio 
waves from a given source arrive at the different receivers of 
an interferometer having traversed different paths through 
the ionosphere, acquiring different phase delays if the density 
is non-uniform. An initially planar set of wavefronts thus be¬ 
comes distorted. This can cause the source to appear shifted 
from its true posi tion, its pea k amplitude to vary and/or its 
shap e to distort (ISmithlll952l : ISooelstr'^l 19971 : lintema et al.l 
I 2 OO 9 II . Multi-path propagation through a highly irregular 
plasma can lead to self-interference of coherent phase fronts, 
causing a dif fraction pattern to develop and giving rise to 
scintillation dMever-Verned [l980l l . Phase delays associated 
with water vapour in the troposphere, although significant 
at high (>1 GHz) frequencies, are ~100 times smaller than 
those of the ionosphere at the low frequencies of the MWA 
dThompson et al]l200ll '). Furthermore, the aridity of the site 
ensures that tropospheric activity is relatively benign. The 
tropospheric contribution can therefore be largely ignored, 
and we consider only ionospheric effects here. 


2.1 Angular position 

If the TEC along the lines of sight from all receivers to a 
given source is the same, then the only effect (ignoring those 
related to polarisation and frequency dispersion, which we 
do not consider here) is to add a constant phase to all an¬ 
tennas. A constant absolute phase is not measurable by cor¬ 
relation interferometry, and so the recorded visibilities are 
identical to what they would be if there were no ionosphere. 
An angularly uniform TEC screen therefore has no effect on 
radio observations; distortions only arise from density inho¬ 
mogeneities. 

For radio waves propagating well above the local elec¬ 
tron plasma frequency cup = (e^ne/eouie)^^^, the refractive 
index n is given by n « 1 — tjp/2ti;^, where uj ^ ujp \s the 
angular frequency of the radio waves, eo is the vacuum per¬ 
mittivity, and e and me are the electron charge and mass. 
Consider the lines of sight and I 2 from two separate in¬ 
terferometer elements to a given radio source, where the 
elements are spaced by a distance D. The optical path dif¬ 
ference between £i and £2 is 

OPD= f n{x)dx- [ n(x)dx = -ATEC , 

Je, ^ ' Ji2 2ij^eom, 

( 1 ) 


where ATEC = TEC(4)-TEC(£ 2 ). If ATEC / 0, then the 
wavefront appears to be tilted with respect to the original 
and this cau s es the apparent positio n of the source to shift 
llSmithiri952l : [Thompson et al.ll200lh . The phase difference 
between the antennas is 

A</ = </(4)-</>(4) = - -ATEC, (2) 

c dTreocme u 


where i/ = is the observing frequency. In the small- 

angle approximation and assuming parallel lines of sight, 
the resulting angular shift in position is 


OPD __1 1 ATEC 

D eorue D 


(3) 


The negative sign in Equation m indicates that sources 
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refract in the direction of decreasing TEC. The depen¬ 
dence implies that for a given ATEC/D, the refractive shift 
is larger at lower frequencies. 

If the collection of sight lines from all antennas to a 
given source passes through a patch of the ionosphere whose 
TEC distribution is well described by a linear ramp (insignif¬ 
icant spatial curvature), then ATEC is linearly proportional 
to D. Equation m then implies that A9 will be the same 
on all baselines. Substituting values for the various physical 
constants, this yields 

Ad « -i^^VuTEC (4) 

as the angular offset of the radio source in an image syn¬ 
thesised using the array, where A6 is in radians, v is in 
Hz and V±TEC (the transverse gradient of the TEC) is in 
electrons m~®. At mid-latitudes where the MWA is sited, 
medium-scale travelling ionospheric disturbances (TIDs) 
driven by atmospheric wa ves are the most famil iar type 
of wavelike perturbation dThompson et all l200lli . These 
have associated VuTEC val ue s of electro ns m“® 

djacobson fc EricksorJ[l992al f9: iDvmond et al.lBoill 'l , which 
translates to A9 ~ lOarcsec at i/ = 154 MHz. 

As long as the length scales of the perturbation are 
much larger than the physical size of the array, irrespec¬ 
tive of the field-of-view (FoV), and V^TEC is roughly con¬ 
stant within an integration time, the perturbed wavefront 
remains approximately planar over the array and distor¬ 
tions are primarily in the form of angular position shifts 
with minimal changes to source shape or amplitude. The 
physical size of the MWA, which has a longest baseline of 
km (comparable to the m ean free path of neutral par- 
ticlej3 in the thermosphere; Ijacchial Il977l l is small com- 
pared to the wave lengths of TIDs, which can be 100-1000 km 
llHunsuckeilll982l ~) . This suggests that neutral gas motions at 
ionospheric heights, including the atmospheric waves that 
drive TIDs, are locally expected to produce only simple 
position shifts of point sources in MWA images. However, 
the energy of these disturbances may potentially drive cas¬ 
cades of sma ller-scale irregula r ities through various plasma 
instabilities llWhiteheadlll97il : lKlostermeverj|l978l : iRottgeil 
S). Plasma turbulence occurring on scales below 3 km 
down to the ion gyrora dius (~5 m in the thermosphere; 

I Woodman fc Ba^]l978l 'l may prod uce higher-ord er wave- 
front perturbations over the MWA (lBookerj|l979l l. leading 
to source broadening and shape distortions. 

2.2 Flux density 

While position shifts are governed primarily by the first 
spatial derivative of the TEC, scintillation and angular 
broadening effects a r e related to higher- order derivatives 
dhee &; Jokipiil Il975l : iMever-Vernetl Il980l 'l. The scattering 
and diffusion of ray paths upon propagation through an in¬ 
homogeneous medium causes focusing and defocusing of ra¬ 
dio waves. Curvature of an initially planar wavefront causes 
otherwise point-like sources to acquire a finite coherence 
length, leading to an increase in apparent width and a de¬ 
crease in peak amplitude. Temporal variations in the sig¬ 
nal may come from intrinsic time variability in the plasma, 

^ Mostly N 2 , O 2 , O and He at the relevant altitudes. 


and/or relative motion between the plasma and the lines 
of sight to celestial sources. This relative drift arises both 
from the rotation of the Earth (~20ms“^ drift of sight 
lines through the ionosphere at the latitude of the MWA, as¬ 
suming an ionospheric altitude of 300 km), and also neutral 
winds causing the bul k plasma to drift at 10-100 ms“^ with 
respect to the ground dPavies et al.l[l973 f. Small systematic 
deviations in flux density may also arise from absorption due 
to collisions between electrons and ions/neutrals, this being 
~1% in the daytime a nd ~0.1% at night for typical MWA 
operating frequencies (iThompson et al. 1200 ill . 

For perturbations to the phase fronts by an amount 
less than ~1 radian within a patch of some critical length 
L, amplitude distortions are small (fractional variations 
much less than unity). This is the regime of weak scintil¬ 
lation, inhabited by the ion osphere most of the time for 
frequencies around 100 MHz llNaravarJll992ll . The timescale 
of amplitude variability is L/v, where v is the relative 
drift velocity between the i onosphere and the sight lines 
dNaravan fc GoodmanI Il989lf . The coherence length L is 
given by the larger of either the baseline length or the Fres¬ 
nel scale rp = %/Xi, the latter of which describes the size of 
the region on the ionosphere from which Huygen wavelets 
arrive at a p oint on the ground approximately in phase 
(ICronvnlll972h . Here A is the observing wavelength and 2 
is the altitude of the ionosphere. Baseline lengths for the 
MWA (several hundred metres) are comparable to or below 
the Fresnel scale at 154MHz (~lkm), implying a charac¬ 
teristic timescale of 10-100 s for amplitude variability due 
to the ionosphere. 

If deflection angles are sufficiently large, rays from a 
given celestial source intersect before they reach the ground, 
producing a fully-developed Fraunhofer diffraction pattern 
that gives rise to large and rapid amplitude fluctuations 
(fractional variations of order unity). This regime is known 
as strong scintillation, and occurs when there are signifi¬ 
cant phase fluctuations on spatial scales below rp. If such 
irregularities are produced by a turbulent cascade driven by 
motions of the neutral atmosphere, then we might also ex¬ 
pect such events to be accompanied by significant position 
fluctuations induced by longer-wavelength modes within the 
inertial range (since outer scales would exceed the neutral 
mean free path of 3km, and thereby the size of the MWA). 
Since rp oc \l\fv and At)) oc 1/iz, when v is small irregu¬ 
larities of a given scale size are both associated with larger 
phase fluctuations and have a greater chance of falling below 
rp. Scintillation events are therefore more likely to occur at 
low frequencies. 

A source whose angular size exceeds Ljz will not scin¬ 
tillate, because the diffraction patterns associated with dif¬ 
ferent parts of the source will blend together and average 
out. For the MWA at 154 MHz, the critical angular diam¬ 
eter below which a source can scintillate, given by the an¬ 
gle subtended by the Fresnel scale at ionospheric heights, is 
about 10 arcmin. Such a source would be resolved in the im¬ 
ages (synthesised beamwidth ~2arcmin at 154MHz). The 
vast majority of radio sources in the sky are not resolved 
by the MWA, and so ionospheric scintillation events, when 
they occur, would be expected to affect a large number of 
sources in the data. This contrasts interplanetary and in¬ 
terstellar scintillation, which are only expected to affect the 
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most compact sources dThompson et al]|200ll : iKaplan et alJ 

I2OI5I) . 

One way to distinguish ionosphere-induced variability 
from intrinsic variability is to check whether neighbouring 
sources are varying with a similar modulation index on a 
similar timescale. However, if it so happens that a highly 
localised patch of irregularities conspires to affect only an 
isolated source (thus mimicking intrinsic variability), it is 
unlikely that this patch will remain in front of the source 
for very long. Because of the rotation of the Earth, in-situ 
stationary ionospheric irregularities will drift with respect to 
celestial sight lines at a rate of about 0.2degmin“^. Given 
an average source spacing of 0.5-1° for a 2-min integration 
with the MWA, it will only take a few minutes a given source 
to drift out of its patch of ionosphere and into its neigh¬ 
bour’s. Variability of an isolated source that continues for 
more than several minutes is therefore more likely to be in¬ 
trinsic (or perhaps of interplanetary/interstellar origin) than 
ionospheric. 

Unlike position offsets, which are controlled by the 
largest-scale density structures, amplitude variations are 
governed by irregularities on the smallest scales. Large-scale 
density variations associated with TIDs and whistler ducts, 
the t wo most common types of structures appearing in MWA 
data dLoi et al.ll2015a| L produce only negligible flux density 
ch anges: for t he rela tively extreme TEC gradients measured 
bv iLoi et al.l (l2015bl i for whistler ducts, with angular deflec¬ 
tions of A9 ~1 arcmin over L ~ 10 km transverse scales, the 
associated fractional amplitude variations are estimated to 


AS _ 

~ ~ ~ir 


(5) 


where S is the flux density. Eor TIDs, which have AO ~ 
lOarcsec and L ~ 100 km, we find that AS/S ^ 0.01%. 
This is well below instrument sensitivities for 2-min integra¬ 
tions (~50mJy; this has contributions from thermal noise, 
and also classical and sidelobe confusion) for typical source 
brightnesses (~0.5 Jy) at 154 MHz. 

A more sophisticated treatment is required for small- 
scale irregularities. Theoretical work considering the effects 
of scintillation on interferometric visibilities has been con¬ 
ducted by a number o f authors, und e r various assumptions 
and approximations ( Cron^ 197^ Goodman fc NaravarJ 


1 19891 : iKoopmansI I201C : IVedantham fc Koopmansll2014| ~). A 


thin phase-screen model is often used, and the spectr um of 
small-scale fluctuatio ns assumed to be a power law dRinol 
1 19821 : Iwheeloni [2003h . The results have been shown to de¬ 
pend on many quantities, including the strength, spectral 
index and inner and outer scales for the phase fluctua¬ 
tion spectrum, the channel width, integration time and fre¬ 
quency of the observations, the altitude and thickness of 
the phase screen, and the FoV and geometry of the ar¬ 
ray. Gompared to position offsets, which depend on just two 
free parameters {v and VxTEG), it is far less straightfor¬ 
ward to estimate the variability of peak flux densities in 
th e hnal synthesised images . The recent theoretical work 
of IVedantham fc KoopmansI ll2014h for realistic turbulence 


^ We arrived at Equation lO under simple geometric optics con¬ 
siderations, which are approximately valid under conditions of 
weak scintillation (no interference effects). 


parameters suggests that ionospheric scintillation may be a 
substantial source of noise variance in the amplitudes and 
phases of visibilities recorded by widefield, low-frequency ar¬ 
rays. Here we empirically assess scintillation effects on source 
flux densities for the case of the MWA at 154 MHz. 


3 OBSERVATIONS, DATA REDUCTION & 

ANALYSIS 

Part of the usual calibration process of radio interfero¬ 
metric data is a removal of ionospheric effects. For con¬ 
ventional instruments, which have narrow helds of view 
and operate at higher frequencies, the ionosphere can 
be assumed to be locally uniform within the FoV of 
each antenna. In this situation, a single phase term 
for each antenna is sufhcient for capturing the effects 
of the ionos phere. Calibration algor ithms such as phase 
referencing (|FomaIont_|£j)erley| 199^) and self-calibratio n 
(IPearson fc ReadheadI 19841 : Cornwell fc Fomalond ll999l L 
which solve for a single phase correction term per antenna, 
are effective at removing ionospheric phases from the data. 

In contrast, the wideheld nature of next-generation in¬ 
struments implies that each antenna sees a large patch of the 
ionosp here, and the assumption of local uniformity b reaks 
down llCotton et al.l[2004l : lLonsdalell2005l : llnteni^l2009ll . The 
ionospheric phase is now a direction-dependent function for 
each antenna, and this cannot be removed using conven¬ 
tional calibration algorithms, which neglect direction de¬ 
pendencies. In this study, we analyse data for which only 
time-independent and direction-independent corrections to 
the antenna phases have been applied. The wide FoV of the 
MWA, combined with its compact size, imply that anten¬ 
nas see almost identical phase screens and so differences in 
the ionospheric phase between antennas are expected to be 
small. The calibration process to which the data were sub¬ 
jected is therefore only effective for removing instrumental 
phase errors (e.g. due to slow temperature-induced changes 
in cable lengths), and leaves the ionospheric phase (domi¬ 
nated by rapid time- and direction-dependent rather than 
antenna-dependent variations) largely untouched. 


3.1 Observing Strategy 

The data were recorded over a 30.72 MHz band centred at 
154 MHz. Data collection began in mid-2013 and took place 
one night each month. During an observing run, which typi¬ 
cally lasted for 9-10 hr, the telescope was pointed along the 
meridian at three different declinations, centred at zenith 
and zenith angles of ~30° {5 = —55?0, —26?7, -|-1?6). A se¬ 
ries of 2-min snapshots were obtained at each declination in 
turn, the cycle repeating every 6 min. The drift of the sky 
during each 2-min scan was accounted for during imaging by 
hxing the phase centre for each snapshot to be at a certain 
RA/Dec. We have sub-divided the data by month and dec¬ 
lination band, such that each dataset corresponds to a hxed 
pointing direction (in Az/El) and contains snapshots from 
a single night. The naming convention here is “year month 
declination”, where for example “2013 Sep —55” refers to the 
dataset obtained in September of 2013, with snapshots cen¬ 
tred at 5 = —55?0 on the meridian. 

This study analyses all MWATS datasets currently 
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available that have passed quality control, namely those for 
which post-calibration gain fluctuations lie within accept¬ 
able tolerances. Although they were not selected based on 
any measure of ionospheric, geomagnetic, tropospheric or 
solar activity, we note that all MWATS observations were 
done at night and happened to be conducted under quiet 
geomagnetic conditions {Kp index < 2). Our results should 
therefore be representative of typical quiet nighttime condi¬ 
tions above the observatory. It should be emphasised that 
these results do not apply to other times of the day, since 
the electron content undergoes large diurnal modulations, 
or observatories at other locations, because of the strong 
dependence on geomagnetic latitude. Ionospheric irregular¬ 
ities are much more prevalent at equatorial and high lati¬ 
tudes, whereas mid-latitude sites such as where the MWA 
is situated are known to experience lower levels of activity 
dFeier fc Kellevl[l98^ 1. 

There are a total of 20 datasets accounting for 10 dif¬ 
ferent nights and a total of about 51 hours of observation, 
spread over a period of slightly more than a year between 
2013 and 2014 (observing dates/times and central coordi¬ 
nates are listed in Table [TJ. We make direct use of the stan¬ 
dard survey images with no special reprocessing. Note that 
MWATS is intended for general-purpose astrophysical tran¬ 
sient and variable searches, and is not designed for charac¬ 
terisation of the ionosphere. However, given that the aim 
of our study is to examine the variability that manifests at 
the time and frequency resolution of MWATS (and not, for 
example, to constrain the structure of the ionosphere), we 
consider the data sufficient for this purpose. 

3.2 MWA Data Reduction 

To calibrate each dataset, we observed a bright source with 
well-modelled emission (a calibrator) for two minutes at the 
phase centre of the instrument. We chose the nearest avail¬ 
able calibrator in declination to the observations. A reference 
image, taken from another low-frequency instrument, was 
used as the starting point for calibration. Time-independent, 
frequency-dependent phase and ampltude calibration solu¬ 
tions were derived based on a frequency-adjusted model of 
each calibrator. Table [T] lists the calibrators used for each 
dataset. 

We used a pre-processing algorithm to flag radio¬ 
frequency interference (RFI), average the data and then con¬ 
vert it into CASA0 (version 4.2) measurement set format. 
RFI flagging was ach ieved using the algorithm AOFLAGGER 
dOffringa et al.l [20121 1. The calibration solutions were then 
applied to the visibilities, which wer e imaged and decon - 
volved using the wsclean algorithm (lOffringa et al.ll2014h . 
An image size of 3072 x 3072 pixels was used with a pixel 
size of 45 arcsec. Stokes I images were formed using a Briggs 
weighting of —1, giving a result closer to uniform than nat¬ 
ural weighting. The wsglean algorithm achieves this by 
forming a complex 2x2 Jones matrix / for each image 
pixel. The images were restored using a circular Gaussian of 
width 130 arcsec, which roughly describes the shape of the 
MWA synthesised beam at zenith, and this was elongated 
appropriately off-zenith to account for the foreshortening 

® http://casa.nrao.edu 


of projected instrument baselines. Primary beam correction 
was performed by inverting the beam voltage matrix B and 
computing , where * denot es the conjugate trans¬ 

pose (details in lOffringa et al]l2014h . The resulting Stokes 
I snapshots were then used in subsequent analysis. It is to 
be noted that residual errors exist in the analytical beam 
model; these are discussed further in i)3.5l 

3.3 Source Extraction and Selection 

We extracted sources from the snapshot im ages using the 
source finder Aegean ([Hancock et al.ir2012| j. This fits a 2D 
Gaussian to each source it detects above a certain noise 
threshold. The outputs from Aegean relevant for this analy¬ 
sis are the best-fit angular position, peak flux density, Gaus¬ 
sian shape parameters and their associated errors. Among 
the extracted sources, we excluded those with negative peak 
flux densities, those with position errors exceeding 1 arcmin 
(as quoted by Aegean, and corresponding to the positional 
uncertainty on a source with S/N ratio ~1), and those for 
which there were fewer pixels above 5cr than the number of 
parameters to be fitted. The purpose of these restrictions 
was to remove the poorer-quality fits. 

We then cross-matched the remainder with either the 
NRA O VLA Sky Survey (NVSS) catalogue dCondon et al.l 
Il998ll or the Sydn ey University Mo longlo Sky Survey 
(SUMSS) catalogue (iMauch et al]|2003l . depending on the 
declination range of the dataseiQ. The cross-matching radius 
we used for most datasets was 1.2 arcmin. However, for the 
datasets 2014 Aug —26 and 2014 Aug—55, where extreme po¬ 
sition offsets were detected (of order the synthesised beam 
or larger), we had to increase the cross-matching radius to 
3.6 arcmin to obtain a cross-matching efficiency (fraction of 
Aegean sources with a match in the external catalogue) 
comparable to the other datasets (~95%). 

To reduce the fraction of spurious sources (which may 
be noise spikes or residual sidelobes), we discarded all 
sources that did not have a match in the external cata¬ 
logues. We also placed a restriction on the minimum number 
of snapshots a source had to appear in (on a dataset-by- 
dataset basis) to be considered for subsequent analysis. The 
value of the restriction was 20 snapshots (roughly 75% of the 
maximum number of snapshots that a source at the central 
declination will appear in) for most datasets, but we low¬ 
ered this to 10 for several smaller datasets that had very few 
sources appearing in more than 20 snapshots. This was cho¬ 
sen as a compromise between excluding sources whose true 
positions were less well represented by the time-averaged po¬ 
sitions (see later in 1 13.411 . and retaining a sufficiently large 
sample for statistical analysis. Note that many sources ap¬ 
pear in more than 20 snapshots, particularly for the southern 
MWATS pointing where 30-40% of sources are present in 30 
or more snapshots. The value of the cut on the number of 
snapshots and the resulting numbers of sources analysed are 
stated in Table [2] which also lists the total number of snap¬ 
shots and source occurrences for each dataset. Note that 
this filtering implies that the resulting source sample is not 

NVSS covers & > -40°, while SUMSS covers <5 < -30°. We 
cross-matched datasets having <5 = —26?7,-|-1?6 with NVSS and 
those having 5 = —55?0 with SUMSS. 
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Table 1. Parameters for each dataset. 


Dataset 

UTC start/end 

AWST_^ start/end 

LST start/ 
end (hrs) 

Central RA/Dec 

Central Az/El 

Calibrator 

2013 Sep-55 

2013-09-16 13:30:39 

2013-09-16 

21:30:40 

21.0 

12°, -55° 

180?0, 61?7 

PKS 2356-61 


2013-09-16 21:24:39 

2013-09-17 

05:24:40 

4.9 




2013 Oct-26 

2013-10-16 13:40:39 

2013-10-16 

21:40:40 

23.1 

44°, -27° 

0?0, 90? 0 

3C444 


2013-10-16 21:34:39 

2013-10-17 

05:34:40 

7.1 




2013 Dec-1-1.6 

2013-12-06 13:51:28 

2013-12-06 

21:51:28 

2.7 

87°, -1-2° 

0?0, 61?7 

Hydra A 


2013-12-06 20:09:28 

2013-12-07 

04:09:28 

9.0 




2013 Dec-26 

2013-12-06 13:55:27 

2013-12-06 

21:55:28 

2.7 

88°, -27° 

0?0, 90? 0 

Hydra A 


2013-12-06 20:07:27 

2013-12-07 

04:07:28 

9.0 




2013 Dec-55 

2013-12-06 13:53:28 

2013-12-06 

21:53:28 

2.7 

86°, -55° 

180?0, 61?7 

Pictor A 


2013-12-06 20:05:28 

2013-12-07 

04:05:28 

8.9 




2014 Mar-1-1.6 

2014-03-06 11:20:48 

2014-03-06 

19:20:48 

6.1 

161°, -1-2° 

0?0, 61?7 

Hercules A 


2014-03-06 21:00:48 

2014-03-07 

05:00:48 

15.8 




2014 Mar-26 

2014-03-03 11:32:32 

2014-03-03 

19:32:32 

6.1 

161°, -27° 

0?0, 90? 0 

Hydra A 


2014-03-03 21:12:32 

2014-03-04 

05:12:32 

15.8 




2014 Mar-55 

2014-03-17 11:07:28 

2014-03-17 

19:07:28 

6.6 

143°, -55° 

180?0, 61?7 

Pictor A 


2014-03-17 20:47:28 

2014-03-18 

04:47:28 

16.3 




2014 Apr-1-1.6 

2014-04-28 10:52:00 

2014-04-28 

18:52:00 

9.1 

210°, -1-2° 

0?0, 61?7 

Hercules A 


2014-04-28 20:46:00 

2014-04-29 

04:46:00 

19.0 




2014 Apr-26 

2014-04-28 10:50:00 

2014-04-28 

18:50:00 

9.0 

209°, -27° 

0?0, 90? 0 

Hydra A 


2014-04-28 20:43:52 

2014-04-29 

04:43:52 

19.0 




2014 Apr-55 

2014-04-28 10:48:00 

2014-04-28 

18:48:00 

9.0 

209°, -55° 

180?0, 61?7 

Virgo A 


2014-04-28 20:42:00 

2014-04-29 

04:42:00 

19.0 




2014 Jul-1-1.6 

2014-07-15 10:46:24 

2014-07-15 

18:46:24 

14.1 

285°, -1-2° 

0?0, 61?7 

Hercules A 


2014-07-15 20:40:24 

2014-07-16 

04:40:24 

0.0 




2014 Jul-26 

2014-07-15 10:44:24 

2014-07-15 

18:44:24 

14.1 

285°, -27° 

0?0, 90? 0 

3C444 


2014-07-15 20:38:24 

2014-07-16 

04:38:24 

0.0 




2014 Jul-55 

2014-07-15 10:42:24 

2014-07-15 

18:42:24 

14.0 

284°, -55° 

180?0, 61?7 

3C444 


2014-07-15 20:36:24 

2014-07-16 

04:36:24 

0.0 




2014 Aug-1-1.6 

2014-08-26 11:00:48 

2014-08-26 

19:00:48 

17.1 

330°, -1-2° 

0?0, 61?7 

Hercules A 


2014-08-26 20:54:48 

2014-08-27 

04:54:48 

3.0 




2014 Aug-26 

2014-08-26 10:58:48 

2014-08-26 

18:58:48 

17.1 

330°, -27° 

0?0, 90? 0 

3C444 


2014-08-26 20:52:48 

2014-08-27 

04:52:48 

3.0 




2014 Aug-55 

2014-08-26 10:56:48 

2014-08-26 

18:56:48 

17.1 

329°, -55° 

180?0, 61?7 

PKS 2356-61 


2014-08-26 20:50:48 

2014-08-27 

04:50:48 

3.0 




2014 Oct-1-1.6 

2014-10-07 11:15:12 

2014-10-07 

19:15:12 

20.1 

15°, -1-2° 

0?0, 61?7 

3C444 


2014-10-07 21:09:12 

2014-10-08 

05:09:12 

6.0 




2014 Oct-26 

2014-10-07 11:13:12 

2014-10-07 

19:13:12 

20.1 

15°, -27° 

0?0, 90? 0 

3C444 


2014-10-07 21:07:12 

2014-10-08 

05:07:12 

6.0 




2014 Oct-55 

2014-10-07 11:11:12 

2014-10-07 

19:11:12 

20.0 

14°, -55° 

180?0, 61?7 

Pictor A 


2014-10-07 21:05:12 

2014-10-08 

05:05:12 

6.0 





“ AWST = Australian Western Standard Time (UTC +8) 


complete. Rather, it is intended to be a reliable sample of 
sources that can be used to probe ionospheric fluctuations. 

3.4 Measuring Position Offsets 

We measured position offsets by computing the angular dis¬ 
placement of each source (uniquely identified by its NVSS or 
SUMSS catalogue name) from its time-averaged position in 
a given dataset. This yielded a set of displacement vectors as 
a function of time for each source. Use of the time-averaged 
position as a reference subtracts away fluctuations that are 
static with respect to the celestial sky. Ionospheric fluctu¬ 
ations would tend to be fixed to the terrestrial sky, while 
fluctuations hxed to the celestial sky are more likely to be 
imaging artefacts (e.g. caused by sidelobes of bright sources). 
In the absence of other causes of time-dependent fluctuation, 
any statistically signihcant scatter in the measured positions 
occurring within the time span of an observation can be 


attributed to the ionosphere (tropospheric delays are only 
~1% those of the ionosphere and can be neglected). A sin¬ 
gle snapshot image from one of the MWATS datasets, along 
with the associated distribution of angular offset vectors, 
is shown in Fig. [T] (further inspection of the temporal be¬ 
haviour reveals that a small-scale TID is passing overhead). 
Organised density fluctuations are widely observed in the 
remaining datasets, with TIDs and fiel d-aligned ducts be- 
ing among those visually identifiable (cf. lUoi et ah I l2015a|ia . 
However, it is not our intention to pursue the physics of the 
underlying phenomena here, but to focus on the broad sta¬ 
tistical properties of the fluctuations that are relevant to 
astrophysical observables. 

The largest -scale TIDs, with periods of ~lhr 
llHunsuckeilll982 h have wavelengths larger than the MWA 
FoV and would be removed by conventional calibration ap¬ 
proaches (these are only effective at removing ionospheric 
fluctuations larger than the FoV). The greatest contributors 
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Table 2. Snapshot and source numbers for each dataset. The third column lists the minimum number of snapshots 
a source in a particular dataset was required to appear in to be included for further analysis. The fourth column 
refers to the number of distinct SUMSS or NVSS sources matched to the candidate sources extracted by Aegean, 
subject to their appearing in at least the number of snapshots listed in the third column. The last column refers to 
the total number of appearances of these sources in any snapshot. 


Dataset 

Snapshots 

Snapshot cut 

Sources 

Source occurrences 

2013 Sep -55 

76 

20 

2923 

75514 

2013 Oct-26 

76 

20 

2894 

67293 

2013 Dec-1-1.6 

56 

10 

2780 

38240 

2013 Dec-26 

58 

20 

860 

19851 

2013 Dec-55 

60 

20 

1239 

33800 

2014 Mar-1-1.6 

55 

10 

2431 

28350 

2014 Mar-26 

54 

10 

5069 

63819 

2014 Mar-55 

52 

10 

2142 

37845 

2014 Apr-1-1.6 

94 

20 

1264 

27467 

2014 Apr-26 

94 

20 

2490 

57154 

2014 Apr-55 

96 

20 

550 

13858 

2014 Jul-1-1.6 

96 

20 

785 

17046 

2014 Jul-26 

96 

20 

2062 

47742 

2014 Jul-55 

96 

20 

1806 

49244 

2014 Aug-1-1.6 

96 

20 

1343 

29588 

2014 Aug-26 

96 

20 

3203 

75593 

2014 Aug-55 

96 

20 

3466 

98358 

2014 Oct-1-1.6 

96 

20 

943 

20886 

2014 Oct-26 

96 

20 

1719 

40221 

2014 Oct-55 

96 

20 

2049 

57397 


to any residual position offsets in MWA data (following stan¬ 
dard calibration approaches) would be those with spatial 
scales between 10-100 km. These include small- to medium- 
scale TIDs, and whistler ducts. The associated timescales of 
position fluctuation due to these irregularities are of the or¬ 
der minutes to several tens of minutes, significantly shorter 
than the restriction on the observing duration of a source 
for analysis (20 snapshots = 2hr). Measuring offsets with 
respect to the time-averaged position thereby isolates these 
short-term fluctuations. A second reason for choosing to use 
the time-averaged position rather than the catalogue posi¬ 
tion as the reference is to cancel out the effects of calibration 
errors in MWATS d ata, which appear as global offsets (dis¬ 
cussed previously bv iLoi et al.ll2015al '). 

3.5 Measuring Flux Density Variations 

As argued in i)2.2l flux density variations associated with 
ionospheric scintillation are expected to occur on short 
timescales, generally below the snapshot cadence of MWATS 
(6min). Systematic variations of the peak amplitudes on 
much longer timescales are likely to be either intrinsic to 
the source, or a result of instrumental/imaging effects. In 
all datasets, we detected strong systematic variations in the 
peak flux densities of unresolved sources as they drifted 
through the FoV, the fractional attenuation being largest 
near the edges and smallest near the centre, forming a 
concave-down modulation pattern that was static with time 
and present for all datasets. Figure[2^ illustrates this modu¬ 
lation pattern. This is likely to reflect a residual error in 
the model of the primary beam used to correct the im¬ 
ages. However, we discuss possible explanations involving 
propagation-related effects later in il5.2l 

To remove the concave-down modulation pattern, we 
applied a high-pass filter to the radio light curves of each 
individual source. We did this by subtracting a smoothed 


version of each light curve from itself, where the smoothing 
interval was set to 9 snapshots (54 min). The ensemble statis¬ 
tics for raw and high-pass filtered light curves are shown in 
Fig. [2] for all sources seen one of the datasets. The high- 
pass filtering operation isolates the short-term variations, 
which we analyse in 1 14.21 Note that Fig. [2] does not directly 
show the flux density S as a function of time, but rather 
the fractional variation in S with respect to the mean value 
of each source, as a function of time expressed in terms of 
the hour angle. Given that the pointing is fixed along the 
meridian, this equivalently measures E-W position across 
the FoV. Light curves for different sources are overlaid on 
top of one another through this transformation. 


4 RESULTS 

In this section, we present our findings on the behaviour of 
angular position offsets and amplitude variations of the sam¬ 
ple of sources selected for this study. We did not investigate 
broadening or shape distortion effects in detail; suffice to say 
that we did not observe severe distortion or smearing effects 
on any of the datasets. The shape parameters reported by 
Aegean for the 2D Gaussian fits to each source were very 
consistent within and between nights, with axial ratios being 
1.1-1.2 on average. We noted the vast majority (~96%) of 
minor axis fits to lie within 10% of the median value taken 
over each dataset, with fitting errors likely to account for 
the few per cent level of scatter. Major axis fits exhibited 
somewhat greater scatter: ~85% of fits were within 10% of 
the median value for each dataset. This can be explained by 
the systematic variation of the point-spread function (PSF) 
major axis over the FoV (caused by the foreshortening of 
baselines), while the PSF minor axis is expected to be al¬ 
most constant. The number of sources extracted by Aegean 
decreased systematically with increasing noise level, consis- 
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UTC 2014-04-26 15:56:00 



CC 

O 

m 


Right Ascension (J2000) 


UTC 2014-04-28 15:56:00 



Figure 1. Top panel: An example 2-min snapshot image from 
2014 Apr —26, zoomed in to the central quarter so that individual 
sources are visible. The full image contains of order 10^ sources 
above 5 (t. Bottom panel: The vector field of angular displacements 
for that snapshot, shown over the whole FoV for sources appear¬ 
ing in at least 20 snapshots. The middle of each arrow marks the 
location of a point source, and the arrow represents the displace¬ 
ment from its time-averaged position, with arrow lengths scaled 
to 150 times the actual displacement distance. The x-ax.is points 
west, the y-axis points north, and (0, 0) marks the location of the 
zenith (a S anson-Flamsteed projection scheme has been used; for 
details, see iLoi et al.|[^15all . Blue and red denote arrows with 
positive and negative a:-components, respectively, to aid visual¬ 
isation. Organised patches of motion are clearly evident, corre¬ 
sponding to density structures with length scales of ~50 km for 
an assumed altitude of 300 km. 
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Figure 2. (a) The ensemble statistics of radio light curves (frac¬ 
tional variation from the mean amplitude of the source) from the 
2013 Dec —26 dataset for all sources brighter than 2 Jy and ap¬ 
pearing in at least 20 snapshots, plotted as a function of hour an¬ 
gle. The pointing is fixed at zenith, and so this equivalently shows 
the fractional variation in amplitude across the primary beam. 
Roughly 10% systematic variations are observed between the cen¬ 
tre and FWHM of the beam. Very similar attenuation patterns 
are seen for all datasets, (b) The resulting light-curve statistics af¬ 
ter applying a high-pass filter (removing >1 hr trends), plotted on 
identical axes. The high-pass filtered light curves are those used 
in subsequent analyses. In both panels, points are coloured by the 
local density on the page to show the shape of the distribution. 


tent with what one might expect from simple considerations 
(fewer sources can be seen when noise levels are higher). 
This suggests that broadening/distortion effects due to the 
ionosphere are insignificant in these datasets and have little 
effect on the performance of the source finder. 


4.1 Position Offsets 

The angular position offsets display a linear dependence on 
A^, as shown in Fig. [3] for the 2013Dec—26 dataset. Re¬ 
fractive behaviour is expected to follow a A^ proportionality 
(see Equation 0, but here we see that the intercept is non¬ 
zero. This may be explained by a ~3 arcsec position fitting 
error, associated with a combination of thermal/confusion 
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Figure 3. The magnitude of angular position offset of sources in 
four frequency sub-bands of the 2013 Dec —26 dataset, plotted as 
a function of squared wavelength. Only sources appearing in at 
least 20 snapshots have been included in the analysis. The line is 
a linear fit to the four points. The error bars on each data point 
associated with random fitting errors are smaller than the marker 
size. 

noise and inaccuracies in the reference position, that is un¬ 
correlated with the ionospheric errors. In addition, multiple 
factors with different A-dependencies may be contributing 
to the observed offsets, e.g. the wavelength dependence of 
the size of the PSF (oc A) and the sky noise temperature 
(oc A^ ®), which affect the amplitude of scatter in position 
measurements. However, our data do not have sufficient frac¬ 
tional bandwidth to decompose out possible contributing 
factors or rule out alternative functional dependences. 

Let us denote the angular position offset of a source 
(measured with respect to the time-averaged position) by 
the vector s, and let its magnitude be s = |s|. The s-values 
measured for two representative datasets are shown in Fig. [I] 
The dashed line marks the expected magnitude of offset if 
the measured displacements were a consequence of Gaussian 
fitting errors. We estimated the fitting error for a sou rce of 
a certain signal-to-noise ratio SNR by llCondonlll997l 'l 


SNRV8hr2 ’ 

where Of, is the width of the synthesised beam (ISOarcsec 
at 154MHz), and we took SNR to be the ratio between the 
peak flux density of the source and the local root-mean- 
square (RMS) noise. 

The observed offsets in all datasets are systematically 
larger than what one might expect from fitting errors alone. 
They typically lie in the range 10-20 arcsec, consistent with 
the TEC gradients expected for fluctuations driven by at¬ 
mospheric waves (see ini). Offsets seldom exceed the pixel 
size of 45 arcsec, being sub-pixel around 99% of the time 
in most datasets. For the July and August 2014 datasets, 
higher levels of ionospheric activity have led to a reduction 
of this rate to 98% in July and as low as 89% in August. 
Mean values of s and the fraction of sub-pixel offsets for 
each dataset are listed in the second and third columns of 
Tabled 



1/SNR 


Figure 4. Magnitudes of the position offsets for sources in the 
(a) 2013 Dec—26 and (b) 2014 Oct-1-1.6 datasets, measured with 
respect to the mean position of the source, against the reciprocal 
of the S/N ratio. The dotted line indicates the expected offset if 
these were due to fitting errors. Points are coloured according to 
the local density on the page to show the shape of the distribution. 
Only one in five points is plotted. The inset figures show the 
histograms of the distributions projected onto the vertical axis. 


Scatter plots of the position offsets are shown in Fig. [S] 
for three representative datasets. The axes of the plot do 
not correspond to offsets in RA and Dec but rather offsets 
along geographic E-W and N-S, computed as described in 
iLoi et aP ( 2015ah . These distributions for the vast majority 
of datasets show noticeable anisotropy, implying a preferred 
direction of source displacement. They tend to be elongated 
into the first and third quadrants, suggesting that fluctua¬ 
tions are preferentially along geographic NW-SE. Elongation 
of the MWA synthesised beam is insufficient to account for 
this bias, a point that we discuss in more detail in 1 )5.11 If 
ionospheric, this indicates that on a long-term basis, there 
is a statistical preference for TEC gradients to be steeper 
in the NW-SE direction. However, at any one time, MWA 
data can exhibit patches of fluctuation oriented in arbitrary 
directions (e.g. the waveform shown in Fig.[T}r displays some 
patches of motion directed NE-SW). 

The anisotropy in the position offset vectors can be 
quantified by observing how the scalar projection se of the 


© 0000 RAS, MNRAS 000, 000-000 




















































Ionospheric effects on time-domain astrophysics 11 


(a) 2013 Dec-26 




(C) 2014Jul+1.6 



Figure 5. Position offsets for sources in the (a) 2013Dec—26, 
(b) 2014 Mar+1.6 and (c) 2014 Jul+1.6 datasets, measured with 
respect to the time-averaged position of the source. The x axis 
points west and the y axis points north. The dotted lines mark 
the location of the origin. Points are coloured according to the 
local density on the page to show the shape of the distribution. 
Only one in five points is plotted. 

offset vectors changes as a function of the projection di¬ 
rection 0. The mathematical de tails of the steps we took 
to compute se are described in iLoi et all (l2015al ~). Briefly, 
we projected the displacement vectors s onto a given direc¬ 
tion 0 on the sky by taking the inner product with a unit 
vector pointing in that direction, to obtain a scalar ampli¬ 
tude Se- The variation of this amplitude (averaged over all 



Figure 6. Position fluctuation power (squared projected offset) 
as a function of the angle of projection 0, measured east of north, 
for the three datasets shown in Fig. [5] The power has been nor¬ 
malised to the number of data points (number of occurrences of 
any source in any snapshot) in the respective dataset. 


source occurrences) with direction is shown in Fig. [5] for 
the three datasets of Fig. [5] In analogy with the electric 
field vectors of light being passed through a polarising filter, 
the quantity plotted is the average squared amplitude (sg) 
(cf. intensity). The preferred fluctuation direction is given 
by the angle 6 ^max at which (sg) maximises, and the degree 
of anisotropy (cf. linear polarisation) is given by 

. Tf — TTl / 2, . / 2, 

^ " M + m ’ ^ ^ max(sg) , m = mm(s 0 ) (7) 

for a given dataset. The values of A and 0max for each dataset 
are listed in the fourth and fifth columns of Table |3l 

Figure [7] which plots the dependence of A on 0max, 
shows that all but one of the 20 datasets have 0 max in the 
NW-SE quadrants. The single dataset with 0max in the NE- 
SW quadrants has the lowest anisotropy (3%), meaning that 
0max is the least well defined. Discounting this one dataset, 
there is a clear preference among all the remaining datasets 
for position offset vectors to be aligned NW-SE (—90° ^ 
Smax < 0°). Furthermore, there appears to be a dependence 
of A on 0max, where for the datasets exhibiting the greatest 
anisotropy, the direction of the anisotropy is more strongly 
E-W. The correlation coefficient of the points in Fig. 0 is 
= 0.37, but rises to 0.56 when the outlier at 0max = —90° 
is excluded. 

The characteristic displacement associated with the 
isotropic (cf. unpolarised) component of the position offset 
vector field can be computed as Siso = -^/m (equivalently, 
the standard deviation of the distribution projected onto its 
minor axis). If the ionosphere is the only source of system¬ 
atic anisotropies in the position offset vectors, then Siso is an 
upper bound on the contribution of non-ionospheric sources 
of astrometric error (e.g. fitting errors). This is listed for 
each dataset in the sixth column of Table [S] Observed val¬ 
ues of Siso are around 7-10 arcsec, several times larger than 
the expected size of fitting errors for typical image sensitiv¬ 
ities and source brightnesses (~3 arcsec, assuming Gaussian 
noise statistics). 
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Figure 7. The degree of anisotropy of the position offset vectors, 
plotted against the direction of greatest position scatter, for each 
of the 20 datasets. Different symbols have been used to distinguish 
the different declination bands. The vertical dashed line separates 
the NW-SE quadrants (—90° to 0°) from the NE-SW quadrants 
(0° to 90°). 

4.2 Flux Density Variations 

For each dataset, we placed an upper limit on the contri¬ 
bution of ionospheric scintillation by analysing the distribu¬ 
tions of AS values versus S. Here S is the time-averaged 
flux density of a source and AS is the difference between 
a flux density measurement and S (after applying a high- 
pass filter; see ^3.51) . Scatter plots of AS versus S for the 
two datasets exhibiting the highest and lowest amplitudes 
of position scatter are shown in Fig. [8] where both quantities 
are expressed in as multiples of the local root-mean-square 
(RMS) noise. The RMS values were computed by Aegean 
based on the interquartile range over an area 20x20 beams 
in size centred about each source. Each individual AS and 
RMS value is a function of source and snapshot. 

In the analysis described below, all AS measurements 
are treated as independent, without consideration of their 
time ordering. This is in line with our goal here to char¬ 
acterise not celestial but possible terrestrial-based activity 
occurring on timescales comparable to or shorter than the 
integration time of each snapshot. We do not perform a light 
curve-based analysis, which is more relevant to searches for 
celestial-based variability (e.g. iGaensler fc HunsteadI [iood : 
iBell et ^|2014^ . It is important to appreciate that the ro¬ 
tation of the Earth causes celestial sight lines to drift with 
respect to the ionosphere at a significant speed (about 1° 
between MWATS snapshots at the same Az/El, larger than 
at 300km altitude), and so the radiation from a given 
source passes through a completely different patch of the 
ionosphere between adjacent snapshots. We therefore do not 
expect scintillation to produce substantial temporal corre¬ 
lations in the light curves of individual sources measured at 
the cadence of MWATS. 

Our approach to separating the scintillation component 
from other contributors to temporal variation relies on the 
assumption that scintillation effects are described by a char¬ 
acteristic fractional variation in amplitude (modulation in¬ 
dex), whereas thermal noise and sidelobe confusion (two im¬ 


portant sources of noise in MWA data) are associated with 
characteristic absolute variations in amplitude. The faintest 
sources will be dominated by absolute variation effects, while 
the brightest sources will be dominated by fractional varia¬ 
tion effects. Assuming that the different types of error are 
independent, the vertical spread as a function of signal-to- 
noise (S/N) ratio is therefore expected to be the quadrature 
sum of a constant and a linear component. 

To decompose these two components, we binned the 
data by S/N ratio into 19 evenly-spaced intervals up to a 
maximum S/N ratio of 100 in steps of 5, with the excep¬ 
tion of the first bin which included all points with S/N ra¬ 
tio less than 10. We computed the variance V of A5/RMS 
values in each bin and then performed a least-squares fit 
to V = V{X), where X was taken to be the mean S/N 
ratio of points in a bin. We assumed the functional form 
V{X) = + n^X^, with and unconstrained pa¬ 

rameters. Absolute-fluctuation effects are quantified by k 
(expressed as a multiple of the local RMS), and fractional- 
variation effects by (which can be interpreted as a modu¬ 
lation index). In all datasets we obtained positive values for 
and implying real values for k and /r and non-zero 
amounts of each type of noise. Values of k ranged between 
0.7 and 0.9, consistent with there being a significant level 
of classical confusion in the images (much of the pixel-to- 
pixel RMS is static), while values of /r ranged between 1-3% 
(these are listed in Table [3] for each dataset). The values of 
/i obtained are an upper bound on the modulation index 
associated with scintillation. We discuss the implications of 
this in > 15.21 We did not find any strong correlations between 
and (s), the median RMS, the median S/N ratio or Smax. 

For this analysis we have chosen to ignore the 
static, concave-down modulation of the peak flux density 
(e.g. Fig. Hi), removing this component from the light curves 
using a high-pass filter. We do not believe that this is due 
to the ionosphere, but rather that it is the residual of an 
imperfect primary beam model. Our quantitative arguments 
against this being ionospheric in origin are presented in > 15.21 


5 DISCUSSION 

5.1 Interpreting Position Offsets 

A possible contribution to the anisotropy in the position 
offset vector field comes from the u, u-coverage of the ob¬ 
servation, which determines the shape of the MWA PSF 
characterising the direction-dependent resolution of the in¬ 
strument. Asymmetries in the PSF can arise from non¬ 
uniformities in the u, u-coverage. If the u, u-coverage pro¬ 
duces a PSF that is elongated along a certain direction, then 
this can increase the scatter in position measurements along 
that direction independently of ionospheric effects. 

We checked to see if the preferred direction of position 
scatter could be accounted for by the PSF shape and as¬ 
sociated fitting errors. We did this by examining the shape 
parameters fitted to each source by Aegean, and compared 
the characteristic position angle for the source fits in each 
dataset with the direction of greatest position scatter. The 
restoring beam is forced to be a circular Gaussian, and so the 
restoration step should introduce no additional anisotropy in 
the shapes of sources. However, sources in general will not be 
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Table 3. Results for the 20 MWATS datasets. The second to last columns, from left to right, correspond to the 
mean displacement magnitude (s) averaged over all sources in the dataset, the percentage /aubpix of displacements 
that are sub-pixel (pixel size 45arcsec), the degree of anisotropy A, the direction of greatest position scatter fimax 
(measured east of north), the amplitude of displacement Sigo associated with the isotropic component of the position 
offset vector field, the upper bound fi on the modulation index associated with scintillation, and the characteristic 
position angle of the PSF (east of north). 


Dataset 

(s) (arcsec) 

/subpix(%) 

A 

0max (°) 

Siso (arcsec) 

h (%) 

bpa (°) 

2013 Sep-55 

9.8 

99.98 

0.18 

-20 

7.3 

1.9 

81 

2013 Oct-26 

13.5 

99.03 

0.44 

-70 

8.7 

2.1 

160 

2013 Dec-1-1.6 

12.8 

99.44 

0.27 

-70 

9.2 

1.5 

113 

2013 Dec-26 

12.0 

99.79 

0.32 

-70 

8.3 

1.5 

38 

2013 Dec-55 

13.6 

99.75 

0.12 

-60 

10.4 

1.4 

85 

2014 Mar-1-1.6 

10.8 

99.90 

0.51 

-55 

6.4 

1.7 

41 

2014 Mar-26 

11.5 

99.95 

0.52 

-70 

6.6 

1.7 

20 

2014 Mar-55 

13.6 

99.84 

0.03 

+15 

10.8 

2.6 

85 

2014 Apr-1-1.6 

10.3 

99.80 

0.19 

-25 

8.0 

1.5 

23 

2014 Apr-26 

9.9 

99.89 

0.26 

-50 

7.3 

1.4 

20 

2014 Apr-55 

11.8 

99.76 

0.41 

-40 

7.6 

1.9 

59 

2014 Jul-1-1.6 

16.4 

97.86 

0.34 

-45 

11.3 

3.4 

149 

2014Jul-26 

15.2 

97.86 

0.48 

-75 

9.4 

2.4 

20 

2014 Jul-55 

16.6 

98.15 

0.41 

-75 

10.5 

2.3 

38 

2014 Aug-1-1.6 

14.7 

97.90 

0.33 

-70 

10.6 

3.0 

160 

2014 Aug-26 

17.2 

94.21 

0.72 

-85 

00 

bo 

2.0 

5 

2014 Aug-55 

22.8 

89.4 

0.64 

-80 

12.5 

2.1 

45 

2014 Oct+1.6 

11.5 

99.86 

0.21 

-30 

8.6 

1.9 

124 

2014 Oct -26 

13.5 

99.62 

0.14 

-90 

10.4 

2.0 

142 

2014 Oct-55 

19.3 

97.72 

0.17 

-40 

14.2 

2.9 

117 


completely deconvolved and so residnal components of the 
flux density will cause sources to be slightly anisotropic, re¬ 
flecting the shape of the underlying PSF. For most datasets 
we observed the distribution of position angles to peak at 
a certain value, which we took to be the characteristic PSF 
position angle for that dataset. These values are listed in the 
last column of Table [S] 

Figure H] plots the direction of greatest scatter against 
the characteristic position angle of source shape fits. It can 
be seen that position angles differ substantially between 
datasets, spanning the full range of values from —90° to 
-1-90°. They are not correlated with the directions of great¬ 
est scatter, which are largely confined between —90° and 
0°. Gaussian fitting errors arising from an asymmetric PSF 
therefore cannot account for the preferential NW-SE offsets, 
which require a separate explanation. 

The observation that the PSF position angles for off- 
zenith pointings are closer to 0° (N-S) than for zenith point¬ 
ings is consistent with the foreshortening of the projected 
baselines when the telescope is pointed towards the north 
or south. This shrinks the u, u-coverage and elongates the 
PSF in the N-S direction. In the case of the zenith point¬ 
ing, it can be seen that position angles tend to be closer 
to E-W. This can be explained by the fact that the array 
is slightly more ext ended in the N-S than E-W direction 
(|Tingav et al.l Eol.ll ). Night-to-night variations in the PSF 
shape are determined by which tiles are flagged/unavailable 
at the time of observation, and also ionospheric conditions. 
This may account for the spread in characteristic position 
angles between observations at the same Az/El. 

The tendency for position offsets to be preferentially 
NW-SE is consistent with reports of similar anisotropies in 
the literature, as detected in airglow measurements of mid¬ 
latitude ionospheric density structures dMartinis et al.|[2006l : 
IPimenta et al.ll2008l ). It has been noted that density struc¬ 


tures in the southern hemisphere tend to be preferentially 
elongated NE-SW, giving rise to TEC gradients steepest 
along NW-SE. This agrees with the statistical preference 
for NW-SE offsets seen in our MWA data. It is speculated 
that an ele ctrodynamic i nstability known as the Perkins 
instability dPerkinj 1 19731 : iMiller etliD Il997l : |h amz ^ Il999l : 
lYokovama et al.l l2^09^ may account for the formation of 
structures elongated preferentially along a certain direction. 
Under this mechanism, vertical undulations (e.g. seeded by 
atmospheric waves) alter the Pedersen conductivity and sup¬ 
porting forces in a self-reinforcing manner, causing an ini¬ 
tial perturbation to be amplified. The growth rate of this 
instability is direction-dependent with respect to the ambi¬ 
ent electric and magnetic fields. However, an investigation 
into whether or not the perturbations detected by the MWA 
indeed form as a result of the Perkins instability is beyond 
the scope of this work. 

That the position offsets match the known preferred di¬ 
rection of ionospheric structures at mid-latitudes, but show 
no correlation with the PSF position angle, supports the 
idea that a significant amount of the scatter in angular po¬ 
sition measurements is indeed due to ionospheric refraction. 
A further 20-30% of the overall scatter may be explained by 
Gaussian fitting errors. 


5.2 Interpreting Flux Density Variations 

In quantifying extrinsic amplitude variability, we have re¬ 
lied on the assumption that the sources used in this study 
are intrinsically non-variable. While this is true of the vast 
majority of celestial radio sources, it is of course possible 
for some to produce transient and/or variable radio emis¬ 
sion. Transient sources, specihcally those appearing for less 
than an hour, would have been excluded by our decision to 
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S/N ratio 



S/N ratio 


Figure 8. Scatter plots of the difference between measured and 
mean flux density AS versus the mean flux density S (each 
axis normalised by the local RMS), for (a) 2013 Sep—55 and (b) 
2014 Aug—55, the two datasets with the lowest and highest aver¬ 
age amplitude of position fluctuations, respectively. The horizon¬ 
tal axis is shown on a logarithmic scale. Points are coloured by 
the local density on the page to show the shape of the distribu¬ 
tion. Only one in five points is plotted. Note that the taper (drop 
in source counts) for S/N ratios below 10, despite faint sources 
being much more common than bright sources, is due to a sample 
selection effect: we required that sources appear in a large number 
of snapshots, biasing our sample to exclude fainter sources that 
might be lost in regions of high image RMS. 

analyse only sources appearing at least 10 or 20 snapshots. 
Individual variable sources may have been present, but these 
are neither the focus of detection in this study nor likely to 
affect the measured bulk trends of the population, and so 
we have ignored their possible existence. 

Physical conditions for the formation of scintillation- 
inducing irregularities are likely to occur over regional-scale 
patches, and so scintillation events are expected to affect 
large numbers of celestial sources in the FoV at once. This 
is an assumption behind our approach to quantifying iono¬ 
spheric effects by examining bulk trends in variability. How¬ 
ever, we cannot rule out the possibility that scintillation- 
causing irregularities form in localised (~1 km-wide) patches 
of the ionosphere, thereby only affecting isolated sources at 
any one time. We have not inspected our data for evidence of 
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Figure 9. The direction of greatest angular position scatter, plot¬ 
ted against the characteristic position angle of Gaussian fits to 
sources in the dataset, for each of the 20 datasets. Different sym¬ 
bols have been used to distinguish the three declination bands. 
Recall that <5 = —26? 7 corresponds to a zenith pointing, while 
5 = -|-1?6, —55?0 are pointings ~30° off zenith towards the north 
and south. A horizontal line has been drawn at 90° separating 
the NW-SE quadrants (—90° to 0°) from the NE-SW quadrants 
(0° to 90°). 

this. In any case, without extensive follow-up observations, 
it would be difficult to distinguish such cases from intrinsic 
variability. The search for intrinsic variability in MWATS is 
the subject of separate ongoing work. 

The finding that fi lies between 1-3% indicates that 
besides thermal noise and sidelobe confusion, there is an 
additional contribution to the amplitude variation of each 
source of order several per cent the source brightness. Some 
of this could be ionospheric scintillation, but gain calibra¬ 
tion errors and residual primary beam attenuation may also 
enter as fractional variation effects. The /r-values are there¬ 
fore an upper bound on, rather than an unbiased estimate 
of, the contribution from scintillation. Decomposing these 
down into the various possible contributors is a task we de¬ 
fer to future work. For a 1-Jy source, /r ~ 1-3% corresponds 
to a 10-30 mjy additional fluctuation, which is comparable 
to the typical image RMS in MWATS. Most sources are 
fainter than this in which case thermal noise and sidelobe 
confusion effects will be more important, but for brighter 
sources this contribution will be noticeable and may need to 
be appropriately incorporated into analyses of light-curve 
variability. 

The bound on /r enables us to constrain the diffrac¬ 
tive scale Tdiff, which can be thought of as the character¬ 
istic length scale of ionospheric phase fluctuations (physi¬ 
cal distance within which phases vary by 1 rad). We apply 
a number of simplifying assumptions: (1) that the MWA is 
compact enough to assume intensity (i.e. zero-baseline) scin¬ 
tillation^, and (2) that the spectrum of phase fluctuations 
is described by the Kolmogorov power law (iKolmogorovI 

® We justify this by the fact that the MWA is centrally condensed: 
112 out of 128 tiles lie within a 750 m radius of the core, and so 
most baselines and baseline separations are shorter than rp. 
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Figure 10. The predicted fractional deviation in peak flux den¬ 
sity (compared to the true peak flux density) of a source as a 
function of position across the FoV for a zenith drift scan, due 
to a smooth ionosphere located at a constant altitude of 300 km. 
This calculation assumed an integration time of 112 s, an observ¬ 
ing frequency of 154 MHz, a background TEC of lOTECU and 
a PSF width of 130 arcsec, with the source passing directly over¬ 
head. Note that the vertical axis is given in multiples of 10“^. 

A modulation index of much less than unity indi¬ 
cates weak scintillation , in which case we have the result 
that fj, « {rp/rdisf^^ IINaravanlllQQd 'l. As argued in i)2.2l 
the decorrelation timescale for weak ionospheric scintillation 
is upwards of 10 s. The integration time for MWATS is thus 
~10 decorrelation timescales, suppressing /r by a factor of 
up to 3-4. Substituting the values for rp ^ 1 km and the ob¬ 
served ~ 1-3%, we obtain rdiff > 20 km. This is consistent 
with being in the weak scintillation regime {rp < rdiff), and 
also matches the sizes of structures detected in MWA data 
(|Loi et al ] |2015alH . This suggests that the empirical value 
of fj, is physically reasonable. 

Systematic modulations in the amplitudes of celestial 
sources can be produced by propagation effects, but we 
now demonstrate that these are not sufficiently large to ex¬ 
plain the observed concave-down pattern in the light curves 
(Fig. [2^). A still, uniform, constant-altitude ionosphere will 
induce smearing of radio sources as they drift overhead, be¬ 
cause the path length through the ionosphere varies as a 
function of viewing direction. This variation is non-linear, 
implying that VxTEC varies across the FoV. The resulting 
smearing (since the offset is proportional to VxTEC) lowers 
the peak flux density relative to its true value. The estimated 
magnitude of this effect for MWATS observing parameters is 
plotted in Fig. 1101 Although the shape qualitatively matches 
the observed concave-down pattern, the amplitude is about 
three orders of magnitude too small and this effect therefore 
cannot account for observations. 

Consider instead the possibility that the attenuation is 
a result of a propagation effect, such as scattering or ab¬ 
sorption, that has a constant probability per unit length of 
interacting with the incoming signal. Let /(()) denote the 
strength of the transmitted signal relative to the original 
signal for a path at zenith angle C- The light curves indicate 
that /(15°)//(0°) « 0.9, and given that f{() = /(O)®®"^^ we 
can then solve for /(O) = 0.05. This means that for the at¬ 


tenuation pattern to be explained by a constant-probability- 
per-unit-length scattering/absorption effect, the decrease in 
signal amplitude needs to be by a factor of about 20 at 
zenith. This is significantly larger than c an be attributed 
to ab sorption, which is <1% at 100 MHz dThompson et all 
I 2 OOII I . If it were caused by scattering due to small-scale irreg¬ 
ularities, then by conservation of total (integrated) flux den¬ 
sity this should produce highly scatter-broadened sources 
(larger by a factor of 20), which are not observed in the 
data. 

That the concave-down attenuation pattern appears in 
all datasets and persists over entire nights of observation 
makes it unlikely that this is due to extreme ionospheric con¬ 
ditions, of the kind that might cause the observed distortion 
via the above effects. It is more likely that these systematic 
variations are instrumental in nature, arising from imperfect 
knowledge of the MWA primary beam. 

5.3 Practical Implications 

We have established that the typical scatter in angular po¬ 
sition is 10-20 arcsec at 154 MHz. Although dominated by 
ionospheric refraction and several times larger than esti¬ 
mated fitting errors, this is sub-pixel for MWATS the vast 
majority (>99%) of the time and thereby almost always sub¬ 
synthesised beam for the MWA at 154 MHz. This implies 
that ionospheric refraction is unlikely to present a major 
problem for the automated cross-matching of sources for 
time-domain astrophysics with the MWA at this frequency, 
as long as a cross-matching radius of ~l-2arcmin is used. 
However, there were several datasets in which higher levels 
of ionospheric activity causing sources to displace by angular 
distances of order the synthesised beam or more necessitated 
a larger cross-matching radius, which we set to 3.6 arcmin. 

For arrays with longer baselines such as LOFAR or an 
extended MWA, ionospheric effects will be a much bigger 
challenge. Second-order spatial TEC derivatives over the ar¬ 
ray will become significant, leading to refractive shifts in dif¬ 
ferent directions for different baselines. In this case, decorre¬ 
lation and shape distortions could occur on a regular basis. 
Arrays with longer baselines will also be able to resolve a 
greater number density of sources, and cross-matching radii 
as large as 1-2 arcmin will be unrealistic. In such cases, a 
thorough calibration that accounts for the full direction de¬ 
pendence of ionospheric phases may be necessary to achieve 
the objectives of time-domain astrophysics. Fortunately, for 
a compact array such as the MWA, ionospheric effects ap¬ 
pear not to adversely affect the ability to conduct these types 
of studies. 

Our analysis does not include the effects of phase cal¬ 
ibration errors, which manifest as a global shift of sources 
with respect to the underlying coordinate grid. We removed 
these by measuring angular displacements with respect to 
time-averaged positions, which isolates only the short-term 
fluctuations. However, the global shifts induced by calibra¬ 
tion errors can be comparable to those resulting from iono¬ 
spheric refraction, and present an independent source of dif¬ 
ficulty for cross-matching. 

We find that short-term amplitude variations contain a 
component, separate from thermal noise and sidelobe con¬ 
fusion, that is described by a modulation index of /r ~ 1- 
3%. This is an empirical upper bound on ionospheric scin- 
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tillation, and appears to be more important than thermal 
noise/sidelobe confusion for sources brighter than about 1 Jy 
for typical MWA sensitivities. This component may possibly 
arise from effects other than scintillation, but we do not pur¬ 
sue the investigation here. We note that ionospheric activity, 
characterised either by (s) or /r, does not exhibit compelling 
correlations with either the RMS noise or the S/N ratio of 
sources. Thus, at least under the observing conditions of the 
datasets used in this study, the ionosphere does not appre¬ 
ciably affect the point-source sensitivity of the MWA. 

We detected a large-scale, concave-down modulation 
pattern in the peak flux density of sources as they tran¬ 
sited through the FoV. The variation is substantial, of order 
10 %, and is most probably related to an incomplete under¬ 
standing of the primary beam gain pattern rather than an 
ionospheric effect (see arguments in >15.211 . In either case, it 
does not represent intrinsic variability and its ubiquity is 
therefore concerning from the point of view of time-domain 
analyses. Obtaining a better understanding of the charac¬ 
teristics of the MWA primary beam is an area of consid¬ 
erable ongoing work, from both theoretical and emp irical 
viewpoints, d escrib ed in part by ISutinio et al.l (l2014l l and 
iNeben et al.l (l2015l ~). 

Finally, we caution that this is not an exhaustive study 
of all data collected by the MWA. The data analysed in this 
work all happened to be obtained under quiet geomagnetic 
conditions, where the Kp index (quantifying global fluctu¬ 
ations of the Earth’s magnetic field) was 2 or less. Our re¬ 
sults therefore pertain mainly to quiet-time effects of the 
ionosphere. Data have so far not been obtained for MWATS 
under signihcantly disturbed conditions, and so the nature of 
possible worst-case scenarios, including the incidence rate of 
scintillation events that may be associated with storm-time 
activity, remains to be established. 


6 CONCLUSIONS 

We have analysed 20 MWATS datasets encompassing 10 
nights of observations under quiet geomagnetic conditions to 
establish how the ionosphere affects the feasibility of time- 
domain science with the MWA at 154 MHz. The two quan¬ 
tities whose statistical behaviour we examined were angular 
positions and peak flux densities. We found that: 

• Angular positions fluctuate by ~10-20 arcsec, these be¬ 
ing below the synthesised beamwidth >99% of the time; 

• There is a consistent preferred direction for angu¬ 
lar position offsets, this being NW-SE (geographic), which 
matches prior reports of anisotropies observed in ionospheric 
density structures; 

• An upper bound on the modulation index associated 
with ionospheric scintillation is ~l-3%; 

• There is a persistent, concave-down modulation of the 
radio light curves over the MWA FoV, but propagation ef¬ 
fects have difficulty accounting for this. 

It appears that the ionosphere does not adversely af¬ 
fect the feasibility of time-domain science with the MWA at 
154 MHz, at least under conditions similar to those exam¬ 
ined here. A cross-matching radius of 1-2 arcmin is sufficient 
most of the time, but under more severe conditions it may 


be necessary to increase this to 3-4 arcmin. Light-curve er¬ 
ror bars may be made more realistic by taking into account 
the existence of an extrinsic source of amplitude fluctuation 
characterised by a modulation index of several per cent. The 
occurrence rates of extreme events cannot be established 
from our data; this demands a broader investigation under 
a variety of geomagnetic and tropospheric conditions. We 
intend to pursue this in future work. 
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